Two electrons in a strongly coupled double quantum dot: from an artificial helium 

atom to a hydrogen molecule 
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We study the formation of molecular states in a two-electron quantum dot as a function of the 
barrier potential dividing the dot. The increasing barrier potential drives the two electron system 
from an artificial helium atom to an artificial hydrogen molecule. To study this strongly coupled 
regime, we introduce variational wavefunctions which describe accurately two electrons in a single 
dot, and then study their mixing induced by the barrier. The evolution of the singlet-triplet gap 
with the barrier potential and with an external magnetic field is analyzed. 
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I. INTRODUCTION 

There is currently interest in developing means of iso- 
lating spins of individual electrons and coupling them in 
a controlled wa ji24MiSiSiL2iS. This problem is equivalent 
to a formation and controlled dissociation of an artificial 
hydrogen molecule. This dissociation is not achieved by 
the increase of the separation of the hydrogen atoms, but 
by the increase of the tunneling barrier separating the 
dots. In such a process a single two-electron dot breaks 
into two one-electron dots. Hence the analogy to the 
break-up of the helium atom into the hydrogen molecule, 
or to nuclear fission, rather than to a chemical reaction. 
The description of the "artificial fission" process cannot 
be accomplished by the weak coupling approaches, such 
as the Hund-Mulliken method, and requires the treat- 
ment of a strongly coupled electron system. We develop 
such an approach in this work. 

The quantum mechanical study of two electron atoms 
dates back to the work of Born and HeisenbergiS on the 
helium atom. (See 11 for a review on two-electron atoms). 
Similar studies of artificial atoms followed the develop- 
ment of quantum dotai^^^ - nanostructures in which 
the number of electrons can be reduced to a desired value 
(ri e = 0,1,2...) in a controllable manner. In contrast 
to atoms, the confining potential of a quantum dot is, 
to a good approximation, quadratic, so a dot containing 
one electron provides a realization of the exactly soluble 
Fock-Darwin model 12 i 15 i 16 i 17 . 

The problem of two interacting electrons in a parabolic 
potential also admits exact solutions, but only for spe- 
cific values of the oscillator frequency, as was shown by 
TauiiiSiiSiSi. The general case was treated analytically 
by the oscillator representation methodSi and variational 
calculations 2 ^, and studied numerically by the following 
approaches: 'exact' diagonalization using Fock-Darwin 
states 23,24 , integration of the radial motion Schroedinger 
equation after separating the center of mass motion^ 5 , and 



a combination of botbJii. The results were compared with 
experimental data 2 ^ 2 ^ and with the Hartree and Hartree- 
Fock methods 23 . The mean field approaches were applied 
to the two-electron system also in the context of symme- 
try breaking they may induce and subsequent symmetry 
restoration by RPA 27 and projection techniques 28 . This 
theoretical problem, as well as the exact and numerical 
solutions mentioned above, proved relevant for a descrip- 
tion of a Wigner molecule consisting of two electrons in 
a quantum dot 
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containing one electron each 
is equivalent to the problem of an artificial hydrogen 
molecule. A variety of methods were applied here: The 
general case was studied by LSDA 5 , molecular orbital 
calculations 3 ^ and the Hartree-Fock approachiiSSi 3 ^^ 7 . 
refined subsequently by 'exact' diagonalization 3 ^ and 
projection techniques 2 ^ 3 ^^ 7 .. The weakly coupled regime 
was studied analytically by the Heitler-London and 
Hund-Mulliken methods 2 3 . While the analytical results 
by the molecular Heitler-London and Hund-Mulliken 
approaches 3 are very useful, the weakly coupled regime 
does not quantitatively describe the experimental situ- 
ation. In molecular description the starting point are 
two well separated quantum dots. Then, as the distance 
between them is reduced, the electrons start tunneling 
from one dot to another in analogy to a chemical bond 
formation. However, in an actual experiment the double 
lateral dot is defined electrostatically by metallic gates 
located above the two dimensional electron gas^t Here 
the distance between the dots is held fixed and the cou- 
pling between them is controlled by means of the inter- 
dot barrier. When the barrier is zero the electrons move 
freely and our system is a single dot, an artificial he- 
lium atom. When the barrier increases, the single dot 
divides into two, the electrons reconfigure so as to avoid 
the barrier and an artificial hydrogen molecule forms. It 
is difficult to find an analogue of the above procedure in 
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the realm of atomic physics. Nuclear physics, however, 
offers an obvious example - fission of a nucleus. In this 
work we demonstrate that viewing quantum dots as 'arti- 
ficial nuclei' rather than 'artificial atoms' offers also some 
computational advantages. 

Our paper is organized as follows: Section 11 describes 
our model consisting of a two-dimensional parabolic po- 
tential perturbed by a Gaussian barrier running along its 
diameter. In Subsection IIA we briefly discuss the exact 
eigenvectors of the single dot problem, found by Taut^S, 
which arc, however, correct only for specific values of the 
magnetic field, different for each state. In Subsection 
IIB we introduce variational wavefunctions which reduce 
to the exact eigenvectors at these specific magnetic fields. 
We calculate corresponding variational energies and com- 



pare them with exact and numerical values. In Section 
III we switch on the barrier and describe the formation 
of molecular states localized in the two potential minima, 
and the effect of the magnetic field on the singlet-triplet 
gap. 



II. THE MODEL HAMILTONIAN 

Our model Hamiltonian describes two electrons mov- 
ing in the (x, y) plane, confined by a parabolic potential 
with frequency u>o, perturbed by a Gaussian barrier of 
width A and hight Vo, and subject to a perpendicular 
magnetic field B: 



J 



H = -(Vr + L4(rl)) 2 - (V 2 + iA{r 2 )f + ^ 2 (r? + rf) + — ^— + V (e 



e a^I 
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Here A(fi) = (a; c yj/4, -u c Xi/4, 0), lu c = ^ is the 
cyclotron frequency of an electron with effective mass 
to* and charge — e placed in an external magnetic field 
B. The lengths are expressed in effective Bohr radii 
a B = 7p§r (where e is the electric permeability) , whereas 
lj , uj c and Vo in effective rydbergs (lRy = 2 ^f, a i )■ The 
magnetic field points in the negative direction of the z 
axis. 



A. Special exact solutions of the single dot problem 

In this and the following subsection we set Vq = 
and address the single dot problem. Change of variables 



R 



fi+r 2 



?%—f\ separates the relative and center of 



mass motion in the Hamiltonian: H = Hr + H r , where 
Hr is the Fock-Darwin one-particle Hamiltonian. The 
detailed study of the radial Hamiltonian H r is presented 
inSS. Here we only summarize the results which we will 
need in the sequel: The relative motion eigenfunction 
with angular momentum m can be expressed as follows: 



>0) e 



(2) 



where a 
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Aiming at the lowest 



excitations of the radial motion one obtains 




(3) 



(where C m are the normalization constants) and the cor- 
responding eigenenergy: 

i,i m (2+ ItoI) , lS 

(4) 



Em ~ 2 '■ " 1 + 2|to| ■ 
However, the two expressions above are valid only if 

I 



2(1 + 2 m 



(5) 



In particular, for to = we get Eq = ARy provided that 



i, and consequently u c 



The 



corresponding wavefunction reads: 

u Q (p) = C p?(l + V2p)e-£. (6) 

At zero magnetic field this is the ground state wavefunc- 
tion, as it is nodeless. For to = ±1 the energy equals 
E±i = zpHMH + 2 Ry, on condition that a = |. This 

implies that lu c — u> c i = 2 A 
wavefunction follows: 
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The triplet radial 



u 1 (p)=C 1 p'(l + \ -p)e 



(7) 



For example, if lo q — then the triplets have the 

energy 2Ry at zero magnetic field, whereas the lowest 
singlet has the energy 4Ry at ui c = w 3.77 Ry. But 
we are not able to get the exact energies at intermediate 
values of u) c . The subject of the next section will be a 
derivation of accurate upper bounds for these energies. 



B. Variational analysis 

To describe the lowest lying states of radial motion 
with arbitrary angular momentum we suggest variational 



wavefunctions inspired by the form of the eigenfunction 
©■ It does not seem reasonable to change the factors 

pi+\ m \ or e - ^" as it would spoil the behaviour of the 
function at zero or at infinity. The only remaining param- 
eter is the one that multiplies p in the bracket. Therefore, 
we introduce the following family of variational wavefunc- 
tions, labelled by the parameter (3: 



u mtf3 (p) = C mtf3P i + W(l+pp)e- E r. 
The corresponding variational energies follow: 

(3 + c m f3 2 



(8) 



(9) 



where d m = ir(l + |m|), e rn = r(§ + |m|), f m = ir(2 + 
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(2|m| + l)V5 
2d„ 



fc TO = — =■ + 2(\m\ + l)e m , 

(2to 2 +4|m| + 3)d„ 
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The minimum /3 m of this simple function is found to be 
at one of the roots of the quadratic equation: 



(10) 



There arises a question of accuracy of our method. To 
study this problem let us return to the example presented 
in the end of the previous subsection where w = § Ry « 
AmeV in GaAs. Clearly, our variational energies repro- 
duce the exact eigenenergies at the specific values of u> c . 
(In this respect they improve upon the interpolation for- 
mula of Taul^.) To evaluate accuracy at intermediate 
values of lu c we solve the radial eigenvalue problem nu- 
merically, by the Numerov method, and compare the re- 
sults with our variational energies (Fig. The picture 
shows that variational and numerical calculations are in 
very good agreement in a wide range of the magnetic 
field. The m^O energies calculated by both methods 
give almost identical results and are represented by sin- 
gle solid lines. The variational m = energy is even 
more accurate than the numerical result, as it is below 
the numerical value for large magnetic fields. Lower ac- 
curacy of the Numerov method is probably due to the 
fact that m=0 state is the only one which is non-zero in 
the singular point p = 0. 

In Fig. ^ we note the familiar singlet-triplet oscilla- 
tions of the ground stat8i&2£»25i212Sl. They are caused by 
a combination of two mechanisms: First, because of the 
orbital Zeeman term — " c 2 £ ' it is energetically favourable 
for the system to rotate. Nonetheless, it is clear that 
without interaction the ground state would be a singlet 
in any magnetic field. Second, Coulomb interaction acts 
stronger on singlets than on triplets, as in the latter the 




FIG. 1: Comparison of the relative motion energies in a 
parabolic dot with uo = f-Ry calculated variationally and by 
the Numerov method. Solid lines: variational m=0,l,-l,2,- 
2,3,-3 states (increasing order at ui c — O.lRy). Dashed line: 
numerical m=0 state. Other numerical states coincide with 
variational in this figure. 



electrons arc kept apart by the Pauli principle. Thereby, 
the gap between a singlet and a consecutive triplet is re- 
duced below its non-interacting value. Nonetheless, with- 
out the Zeeman term the ground state is always a sin- 
glet no matter how strong is the interaction or the term 
quadratic in We conclude that, in the presence 

of interaction, the increase of the magnetic field results 
in an increase of ground state angular momentum. As 
states with even m are singlets, whereas states with odd 
m are triplets, it causes singlet-triplet oscillations. 



III. DOUBLE DOT - FORMATION OF 
MOLECULAR STATES 



In this section we study the ground state and the first 
excited state energies and densities of two electrons in a 
double dot and their evolution with a magnetic field. For 
this purpose we go back to our model Hamiltonian (|T|) 
and set Vq > 0. Clearly, the potential of the barrier: 



(11) 



couples the motion of the center of mass and the relative 
motion. Nevertheless, we choose variational wavefunc- 
tions as products: 



U m (p, ip, R) = u m {p, Lp)^a{R), 



u m (p,ip) = U m .j3 m (p)- 



'2tt 



(12) 
(13) 
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where the center of mass wavefunction is just the Fock- 
Darwin ground state: 



*o(R) =2^-e~ 2aH . (14) 

The corresponding energy of the center of mass motion 

is E cm — \JojQ + ^f. This choice is justified by the fact 
that the barrier couples only every second center of mass 
wavefunction as a result of parity conservation. 

I 



It is our goal to find matrix elements H m ^ n := 
(U m \H\U n ) . To this end, we calculate the effective po- 
tential V e ff(p,(p) = (Wo|H|^o) which acts only on the 
relative motion coordinates: 



V e ff(p, <p) 



Vl 2 



1 



(15) 



where 7 = 2Ay / a. Now we evaluate the matrix elements 
of the effective potential: 



(u m \V e ff\u n ) = VoC m> j3 m C n ^ n 



(/3 m + /3 n )X( 7 2 + l, 



m 



-,n — m) + PmflnK^ 2 + 1, 



2 

4 + Iml 



, n — m) + 
-,ti- m)), 



(16) 
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FIG. 2: The profile of the double dot potential with uiq 
§ Ry, V = 1 Ry, A = 0.5 a B . 



(l+±cos 2 ip) 



-dip. For some 



where K(x,k,l) = T(k)± 

values of k, I this integral can be expressed in terms of 
elementary and elliptic functions. In practice, however, 
we evaluate it numerically. Finally, we have the desired 
expression for matrix elements of H: 



{E m (p m ) + E cm )5 n 



\V e 



eff\ 



(17) 



We focus attention on the potential with parameters 
ujq = I Ry, Vo = IRy, A = 0.5 a B - (The profile is 
shown in Fig. |2J. 

In the two subsections that follow we study the be- 
haviour of the lowest lying singlet and triplet energies 
and densities at < lo c < 1. The third subsection de- 
scribes the effect of the barrier on the singlet-triplet gap. 



FIG. 3: The lowest lying singlet and triplet energies in a 
double dot with ujq = | Ry, Vo = 1 Ry, A = 0.5 as- Dashed 
lines: Eq - singlet, E±\ - triplet. Solid lines: So, ±2 - singlet, 
-E±i,±3 - triplet. Dotted lines: numerical singlet and triplet 
energies calculated by the configuration-interaction method. 
(In increasing order of energy at u a = 0). 



A. Energies and wavefunctions of singlets 

As a first approximation of the singlet energy we take 
only the uq state under consideration. The resulting en- 
ergy: 



(u |Vef/K) 



(18) 



is plotted as a function of the magnetic field in Fig. [3]and 
compared with numerical calculations performed using 
the configuration-interaction method. 

Our variational energy is about 0.08 Ry higher than the 
numerical one at lu c = and the discrepancy increases 
to 0.15 Ry at u> c = IRy. The reason for that is clear 
from Fig. 2] about lu c = IRy the m = state is close 
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to degeneracy with to = 2, therefore the latter should 
also be taken into account. On the other hand, to = 2 
state is degenerate with m = —2 at lo c — so, in fact, 
we should consider both of them. We diagonalize the 
resulting 3x3 matrix and note that its lowest eigenvalue 
Eq,±2 differs from the numerical result by about 0.05 Ry 
at zero magnetic field (see Fig. |3J) and this discrepancy 
becomes even smaller at larger u> . We recall that m/0 
states vanish when the distance between the electrons 
is zero. Consequently, this increase in accuracy is the 
first manifestation of the formation of molecular states: 
contribution of m = ±2 states pushes the electrons apart 
and, by a non-trivial angular dependence, locates them in 
the dots. To support this statement we calculate the two- 
electron density g in this approximation. Let us denote 
by $(r*i,r2) the two-electron wavefunction. Then: 

g(fi) = 2 J |$(ri,r 2 )| 2 cF 2 . (19) 

We have calculated the wavefunction as a linear com- 
bination of functions U m (p,ip), with the coefficients A m 
found by numerical diagonalization. Recalling that p = 
\fotr, r= (x, y) = ?2 — ?\ we obtain: 

$(n,f 2 ) = y^A m C m>0m — a l ^{x + iy) m ■ 

771 

•(l+/Vo>)e~ Q( ^ +r1) . (20) 

The integral (|19|l is evaluated numerically at cu c = Ry 
and 1 Ry. The results are compared with densities ob- 
tained by the configuration-interaction method. (See Fig. 
0] Fig. [SI. We note that the configuration-interaction 
calculations give a significantly lower density in the cen- 
ter of the dot than our variational method, especially at 
u> c = Q. This can partly be attributed to the fact that 
we neglected the radial motion and center of mass ex- 
citations. The first excited state of radial motion with 
to = is close in energy to the m = ±2 states (see Fig. 
I10|l and, in fact, lowers the density in the center of the 
dot as we show in the Appendix (see Fig. 1 1 lf> . On the 
other hand, the center of mass excitations which couple 
to the lowest singlet are separated by about 0.5 Ry from 
to = ±2 states and we will not discuss them further. 

Both methods predict that the electron density in the 
center of the double dot decreases with the magnetic field. 
Similarly, in the Heitler-London approach^ one observes 
that the overlap between the left and right dot wavefunc- 
tion decreases as a function of ui c . From our perspective 
the explanation of this effect starts from the situation in 
a single dot. There, as the magnetic field increases, the 
lowest singlet has larger and larger angular momentum. 
Consequently, the density has a circular ring-like shape, 
with a minimum in the center, caused by the centrifugal 
barrier. Now, when the inter-dot barrier is switched on, 
the ring shrinks into two peaks as a result of mixing of 
angular momenta and radial excitations. 

Summarizing our discussion in physical terms, the two- 
electron droplet in the lowest singlet state acquires ro- 
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FIG. 4: The two-electron density of the lowest singlet calcu- 
lated by diagonalization using m = 0, ±2 states at lo c = Ry 
(upper fig.), 1 Ry (lower fig.). 
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FIG. 5: The two-electron density of the lowest singlet calcu- 
lated by the configuration-interaction method, at u) c = Ry 
(upper fig.), 1 Ry (lower fig.). 
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FIG. 6: The two-electron density of the lowest triplet calcu- 
lated by diagonalization using m — ±1, ±3 states at uj c = ORy 
(upper fig.), 1 Ry (lower fig.). 
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FIG. 7: The two-electron density of the lowest triplet calcu- 
lated by the configuration-interaction method at uj c — Ry 
(upper fig.), 1 Ry (lower fig.). 



tating and vibrating components, when the barrier is in- 
creased and a magnetic field applied. The vibrations have 
to be included at low magnetic fields to obtain qualita- 
tively correct two-electron density. They counterbalance 
the non-rotating (m=0), peaked in the barrier, compo- 
nent of the droplet. 

B. Energies and wavefunctions of triplets 

Let us now describe the lowest lying triplet. First, 
we remark that the triplet energy, as a function of the 
magnetic field, will have a vanishing slope at uj c = 0. 
In fact, breaking of circular symmetry will eliminate the 
degeneracy of m = ±1 states at w c = 0, making the 
wavefunction $ of the lowest triplet real (up to a con- 
stant complex phase). Linear dependence of triplet en- 
ergy could originate only from the orbital Zccman term. 
But ($\L Z ® I + I ® L z \$) = for any real $. Now 
we illustrate this general argument with a calculation of 
actual triplet energies taking only m = ±1 states under 
consideration. In this approximation the matrix elements 
of the Hamiltonian i|17|l read: 

Hi A = Ei((3i)+E cm + (ui\V eff \ui), (21) 
= E_i[(3i) + E cm + {u-i\V eff \u-i), (22) 
= (ui\V eff \u.i). (23) 

We recall from Subsection IIB that Ei(0) = -^ + F{(3) 

and E-i(p) = !f + F(p), where F((3) = 2a %gffi%j£ . 
They are both minimal for the same value of the varia- 
tional parameter j3 = j3\ which solves the equation l|10|l . 
The eigenvalues of the hcrmitian matrix defined by 1)211 
1231) are readily obtained : 

E± =F + E cm + (u 1 \V eff \u 1 )± 

\/^l + \(ui\V eff \u^W. (24) 

The degeneracy of the states m = ±1 at zero field 
has been eliminated: there is a gap of 2|(-ui|V r e //|u_i)|. 
Moreover, the linear term ±w c /2 is no longer present. In- 
stead, there is a term y^~+ | (t/i | V^// |it_i) | 2 quadratic 

for small lo c . The lowest triplet energy E- = E±\ is plot- 
ted in Fig. 13 At zero field it is only about 0.4 Ry higher 
than the numerical value. As the densities of m = ±1 
states are not peaked in the barrier, it is not a surprise 
that the accuracy of our calculations is better than in the 
previous case. The discrepancy increases, however, with 
the magnetic field. To improve upon our approximation 
we take also m = ±3 states into account. The resulting 
energy -E±i,±3 does not differ much from the previous 
one at uj c — 0, but a remarkable accuracy was achieved 
at higher fields (see Fig. |3J). As before, we plot electronic 
densities at lu c = Ry and 1 Ry. (See Fig. El Fig. [7J. 
The low triplet density in the center of the double dot 
is inherited from the single dot. Loosely speaking, the 



two electron droplet in the lowest triplet state consists, 
from the outset, only of rotating components. After in- 
creasing the barrier and applying a magnetic field it will 
acquire components which rotate faster. As contrasted 
to the lowest singlet state, one can obtain a qualitatively 
valid description neglecting vibrations. This distinction 
between the lowest singlet and triplet states is likely to 
remain valid in more realistic double dots. It may find 
applications in spectroscopic measurements sensitive to 
charge distribution (e.g. in quantum point contact mea- 
surements). 



C. Singlet-triplet transition 

We define the singlet-triplet gap J := E t — E s , where 
E t (E s ) denotes the energy of the lowest lying triplet 
(singlet) and plot the J(w c ) dependence for a single dot 
(ujo = | Ry, Vo = 0) in Fig. |SJ Next, we move to the case 
of a double dot setting the parameters as in the previous 
subsections: u>o = § Ry, Vb = 1 Ry, A = 0.5as and 

choose Et = E±i t ±3, E s — Eo,±2- The respective J(w c ) 
function, plotted in Fig. [51 compares well with numerical 
results which are also presented. 

In Section II we argued that the increase of the mag- 
netic field results in an increase of ground state angular 
momentum leading to singlet-triplet oscillations of two 
electrons in a single dot. Although in a double dot the 
eigenstates do not have a definite angular momentum, 
the orbital Zeeman term and Coulomb interaction are 
still responsible for singlet-triplet transitions. In view of 
this fact and the discussion from the previous subsec- 
tion on mixing of angular momenta, it is not a surprise 
that the barrier only smoothed out the sharp edges of the 
J(w c ) function which were a direct consequence of angu- 
lar momentum conservation. In particular, the positions 
of the first crossing and minimum are similar as without 
the barrier. 

Since the precise shape of the dot does not play a 
role in the above discussion it seems to us possible that 
these qualitative features of the J(lo c ) dependence (the 
crossing, the minimum and the signature of the angular 
momentum conservation) are independent of the confin- 
ing potential. (See e.gml on elliptical dots). This claim 
is also reinforced by a recent analysis by Scarola and 
Das Sarma^ who related the singlet-triplet transitions 
to changes in vorticity of the two-electron wavefunction. 
Although their model potential and variational wavefunc- 
tions differ from ours, their results reflect the features 
mentioned above. Consequently, we expect that also re- 
alistic quantum dots of irregular shape will exhibit these 
properties. On the other hand, more detailed features, 
like the value of J(u> c ), the number of crossings or their 
precise positions will certainly depend on the shape of 
the dot. 
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FIG. 8: The singlet-triplet gap as a function of the magnetic 
field in a single dot with loo = | Ry. Variational calculations. 
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FIG. 9: The singlet-triplet gap as a function of the magnetic 
field in a double quantum dot with loo = \Ry, Vo = lRy, A = 
0.5 as. Solid line: variational results. Dotted line: numerical 
calculations by the configuration-interaction method. 



IV. CONCLUSION 

In this paper we discussed the problem of the transition 
from an artificial helium atom to a hydrogen molecule 
as a function of the barrier potential. We illustrated 
it with a simple calculation of the lowest lying singlet 
and triplet states in the strong coupling limit, where the 
Heitler-London approach is not valid. To achieve this 
goal we introduced variational wavefunctions which de- 
scribe accurately two Coulomb interacting electrons in 
a parabolic quantum dot. The singlet and triplet ener- 
gies, and the singlet-triplet gap J, were calculated as a 
function of the barrier potential and the magnetic field. 
The origin of the singlet-triplet transition was discussed. 
We hope that these variational functions will also be use- 
ful in developing methods of isolating spins of individual 
electrons and coupling them in a controllable manner in 
lateral quantum dots, as well as in other areas of research, 
e.g. to study with analytical expressions the formation 
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FIG. 10: The energy of the first excited state in a single 
parabolic dot with luo = § Ry. Solid line: variational calcu- 
lation. Dotted line: numerical calculations by the Numerov 
method. 



of Wigner moleculesS^S. 

We add two remark of a technical nature: First, in the 
case of a double quantum dot we were using the param- 
eters (3 optimized in a single circular dot. One could as 
well optimize them in the potential under study without 
much additional effort. It turns out, however, that the 
gain in accuracy is very small (at least in the case of the 
double dot presented here) so we did not pursue this ap- 
proach. Second, a more realistic model of a double dot 
would be a single elliptical dot perturbed by a barrier. 
Since an elliptical confinement still separates the center 
of mass and relative motion it is not difficult to adapt 
our method to this case. Quantitative agreement with 
numerical results would require, however, more effort as 
our variational wavefunctions arc optimized in a circular 
dot. On the other hand, our approach works well in the 
case of a single elliptical dot (no barrier) or a quantum 
ring (a single circular dot with a circular barrier). 



the first excited state of the radial motion with m = 0. 
To this end, we proceed as follows: We start from a cer- 
tain exact solution of the single dot problem. Then we 
construct a variational wavefunction of a similar func- 
tional form and orthogonalize it to the ground state vari- 
ational wavefunction. Finally, we diagonalize the double 
dot Hamiltonian using the excited state m — 0* together 
with the previously studied to = 0,to = 2,to = — 2 
states. 

From the analysis by Tau1p2S we obtain that at a — 
there is an eigenstate of energy E = 1 Ry given by: 

u (p) = Cop^(l + Vl2p + 2p 2 )e-£. (25) 

As a matter of fact, at zero magnetic field it is the ground 
state, since it is positive. But the binding potential corre- 
sponding to this eigenstate is much weaker than the one 
we have under study. We expect that at a stronger bind- 
ing potential the first excited state with m — will have 
such functional form so we describe it by a variational 
wavefunction: 

= 6 O ,0P* (1 + h + Sp 2 )e-^, (26) 

where C a is a normalization constant, /3, 5 are varia- 
tional parameters. First of all, we have to make sure 
that the state is orthogonal to Mo,/3 : 

/>oo 

= / % fi(p) u o,p (p)dp = 
Jo 

\{l + ~V^(A> + $) + $ + M + -sfrSpo). (27) 
This implies that S — — Ai — A.20, where 



Acknowledgments 

The authors acknowledge support by the Institute for 
Microstructural Sciences, NRC, by the National Science 
and Engineering Research Council, and by the Canadian 
Institute for Advanced Research. One of the authors 
(W.D.) also acknowledges support from the EC Research 
Training Network 'Quantum Spaces - Non-commutative 
Geometry'. The authors thank M. Pioro-Ladriere, A. S. 
Sachrajda, M. Abolfath, C. Dharma-wardana, and W. 
Jakobiec for stimulating discussions during the course of 
this work. 

APPENDIX: THE EFFECT OF RADIAL 
EXCITATIONS ON THE DENSITY IN A 
DOUBLE DOT 

In this Appendix we calculate the two-electron density 
of the lowest singlet in a double dot taking into account 



_ i + ±0t7? o 

1 1 + fvW 

2 i + fV^/V 



Next, we evaluate the normalization constant: 



(28) 
(29) 





d + 




(30) 


d = 


A? 




(31) 


e = 


1 

2 V 


3 

n - -\/7rAi + 2AiA 2 - A 2 , 


(32) 


f = 


1 

2 " 


- jV^A 2 + Al 


(33) 



Now we are ready to calculate the variational energy: 
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- ~ [°° 11 

E„(J3) = 2a J u oJ (p){ -d* + [-^ + p '- j-j] }u aJ ,(p)dp, 



i = ( _L +v ^ ) _ 2a2(1+ i v /i ) _ ai( ^ + _L )+2iiA2(4+ 3yi; )i 



density 



0.12 r 




FIG. 11: The two-electron density at uj c — i?y calculated 
using m = 0, m = 0*, m = 2, m = — 2 states. 



< fi 0,/§ o l^//|tim> 



(The function K was defined in Section III.) Together 
with matrix elements from Section III we have all input 
necessary to diagonalize the double dot Hamiltonian in 
the subspace spanned by the m = 0, m = 0* , m = 2, m = 
—2 states (where 0* denotes the first excited state with 
m = 0, determined in this Appendix). Having obtained 
the wavefunction, we calculate the density at oj c = 
and plot it in Fig. 1111 Comparison with the density at 
LD C = calculated previously (Fig. indicates that the 



(34) 
(35) 
(36) 
(37) 
(38) 

I 

The optimal variational parameter /3o can be found by 
solving a quadratic equation analogous to l(TU|> . The en- 
ergy is then obtained substituting it back to equation 
(|35|l . At this stage a comparison can be made between 
our variational approach and a numerical solution of the 
single dot problem by the Numerov method. The results 
are plotted in Fig. ^| Similarly as in Section II, we note 
very good agreement between the two methods. Finally, 
we arrive at a task of computing matrix elements of the 
effective potential V e ff. 



(39) 



(40) 

I 

radial motion excitation lowers the density in the center 
of the double dot. 



Vo-fC 2 - 

-_M^(K( 7 2 + 1, 1, 0) + {pi + 2S)K(j 2 + 1, 2, 0) + 
V7 + 1 

~8 2 K{ 1 2 + 1, 3, 0) + 2/3 ^( 7 2 + 1, |,0) + 2P~6K(j 2 + 1, |,0)), 

v olC sC m/9 2+H -, w , , 3 + |m| , 

£P= (K 7 2 + 1, — + (p m + /3 K 7 2 + 1, — =!— 1 , m) 

Vf + 1 2 2 

CflbAn + W7 2 + 1, m) + ^ m K( 7 2 + 1, m)). 
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